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Abstract 

In the early Universe matter was crushed to high densities, in a manner similar 
to that encountered in gravitational collapse to black holes. String theory suggests 
that the large entropy of black holes can be understood in terms of fractional branes 
and antibranes. We assume a similar physics for the matter in the early Universe, 
taking a toroidal compactification and letting branes wrap around the cycles of the 
torus. We find an equation of state pi = Wip, for which the dynamics can be solved 
analytically. For black holes, fractionation can lead to non-local quantum gravity 
effects across length scales of order the horizon radius; similar effects in the early 
Universe might change our understanding of Cosmology in basic ways. 
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1 Introduction 



At the beginning of the Universe we expect that energy density was high, and it is hkely 
that quantum gravity was important. Since string theory provides a consistent theory of 
quantum gravity, we must ask what kind of states are expected in string theory under 
these conditions. Several ideas have been considered, using either string theory or string 
inspired constructions [H [21 [3l H] . 

Another place where matter gets crushed to high densities is in the formation of a 
black hole. In the classical picture of a black hole the curvature is low at the horizon, 
and large at the singularity. If we consider quantum mechanics on such a background 
geometry we run into the the black hole information paradox [5], which implies that 
unitarity of quantum mechanics is lost. 

String theory has made considerable progress in understanding black holes. We can 
understand the entropy of extremal and near extremal holes |6l [3, E], and obtain Hawk- 
ing radiation as a unitary process where excited string states decay [9]. This suggests 
that string theory will change our naive picture of the black hole geometry and allow 
information to leak out in the Hawking radiation. 

Several computations have suggested a 'fuzzball' picture of the black hole interior, 
where the quantum gravity effects are not confined to the vicinity of the singularity, but 
instead spread out all through the interior of the horizon. The key effect is 'fractionation': 
when different kinds of branes are bound together they split up into fractional brane units 
[To] . We can regard the large entropy of the black hole as a consequence of fractionation: 
the entropy calculation just counts these fractional brane units with their appropriate 
spins and fermion zero modes. Fractionation is also responsible for the low energy of 
Hawking radiation quanta. More qualitatively, we can say that the fuzzball picture of 
the black hole interior is also a consequence of fractionation; fractional branes are low 
tension objects that can stretch to horizon scales instead of just planck distance from 
the singularity. The concrete computations leading to the fuzzball picture construct the 
microstates that account for the entropy. For 2-charge extremal holes we can understand 
all microstates, and for 3 and 4 charge extremal cases subfamilies respecting one or more 
U(l) symmetries have been constructed [121 [13] • In each case the microstate is found 
to have been modified in the entire interior of the hole, and there is no horizon. If the 
fuzzball picture were true it would resolve the information paradox, since information 
can escape from the surface of the fuzzball, much like it leaves from the surface of a piece 
of burning coal. 

In this paper we wish to ask the question: can we apply our understanding of black 
holes to say something about the Cosmological singularity? In Fig{T](a) we depict a 
traditional radiation filled Universe. We know that we can get a larger entropy for the 
same energy if we put the mass into black holes of sufficiently large radius; we depict this 
in Fig|T](b). But our Universe does not look like this at all; if black holes had formed 
at early times they would continue to exist till today (unless they were small enough to 
have Hawking evaporated by now). 
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The situation does not change in any material way if we replace the conventional 
picture of the black hole interior with a 'fuzzball' (FiglU^c)); this affects only the interior 
of the hole and not gross properties like the classical attraction between holes. 

But if the maximal entropy state of a black hole is this quantum fuzz, then perhaps 
the maximal entropy state of the Universe is given by such a quantum fuzz filling the 
entire Universe (Fig{T]^d)). Using the microscopic expressions for black hole entropy we 
conjecture an equation of state for this fuzz, and find the evolution of the Universe with 
this equation of state. 





(a) 



(b) 




(c) 



(d) 



Figure 1: (a) Radiation filled Universe (b) All matter in black holes (c) Fuzzball 
picture suggests that interior of horizon is a very quantum domain (d) Quantum fuzz 
filling the entire Universe. 



2 Fractional brane states and entropy 

In this section we will recall some results from the string description of black holes, which 
will motivate our ansatz of the fractional brane state and its entropy. A more detailed 
review of these results can be found in [Ti] . 
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We will work with 10+1 dimensional M theory, using on occasion the language of 9+1 
dimensional string theory when discussing branes. We will let the 10 space directions of 
M-theory be compactified to T^°. We will denote the spacetime dimension as D. 

FigIT](a) depicts a Universe filled with radiation. M-theory has massless quanta, so 
we can certainly achieve such a state. Let us fix the lengths of the sides of the torus, 
and explore the entropy as a function of the total energy. If the spacetime dimension is 
D then 

S ~ (2.1) 
Thus if the Universe was filled with massless radiation we would get 5* ~ for the 
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11 dimensions of M-theory and S ~ Ew for the 10 dimensions of string theory; in the 
latter case has been compactified to a small length so that quanta along x^^ are not 
excited. 



(a) 



(b) 



Figure 2: (a) A string can wind several times around a compact cycle and carry vibrations 
(b) In the 'brane gas' model branes can wrap cycles and carry vibrations. 



2.1 Two charges 

Since we have extended objects in our theory, we can wrap them around the cycles of 
the torus. Consider string theory and let a string be wrapped rti times around a cycle of 
the torus; let the length of this cycle be L. We can add excitations to this string, which 
split up into left movers and right movers. First let the string be a heavy 'background' 
object, with the excitations as small vibrations. The excitations form a massless gas in 
1 + 1 dimensions. The total length of the string is Lt = niL. The energy and momentum 
carried by the left movers is of the form 

El = \Pl\ = -r = (2.2) 
The entropy of the left movers is 

Sl = 2V27T,/Wr^ (2.3) 
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where the dependence S ~ y/riiUp comes from the way the momentum can be distributed 
among different harmonics on the string, and the coefficient arises from the fact that 
there are 8 transverse vibration modes of the string and 8 fermionic superpartners of 
these modes [6j. Note that the Up units of Pl broke up into riinp 'fractional' units of 
momentum because the string itself was a bound state of rii singly wrapped strings; this 
is a simple example of the fractionation mentioned above [TO] . 
Adding in the right movers we have 

S = 2V27T^ {^p + ^p) (2.4) 



Of course we should not really regard the vibrations of the string as small oscillations 
in general, and to carry out the full computation we note that the total energy E oi a. 
string state is given by 

E^ = ih.LT - + SttTNl = {h,LT + + ^tiTNr (2.5) 

ij ij 

where T is the tension of the string, hi = rii — hi, hp = Up — hp give the net winding and 
net momentum carried by the string, and Nl, are the left and right excitation levels. 

The entropy is 

S = 2V27r( + v^) (2.6) 

For vanishing net winding and momentum = 0, hp = we get 

S = 2y^ ^ (2.7) 



This is a faster growth of S than (12. ip . and leads to the well known Hagedorn transition. 

We can understand the above dependence S ^ E also in the more elementary com- 
putation (12.41) . The Universe will have no net string winding, so we will have winding 
as well as anti-winding modes. On the winding mode we have left movers (momen- 
tum) and right movers (anti-momentum), and similarly for the anti-winding modes. 
We find that the entropy is optimized if we put as much energy into string winding 
((^1 + hi)LT = 2niLT = -|) as in the momentum excitations. This gives 

ni = hi E, Up = hp E, ~ y/riiUp ~ E (2.8) 

in agreement with ( 12. 7p . 

The purpose of carrying out the estimate in the crude form (12. 8p is that we wish to 
talk about fractional branes and antibranes. In the count (12. 7p we had a closed loop 
of an excited string, but we see that we can regard the two sides of this loop as string 
'winding' and 'antiwinding', and the excitations as 'momentum' and 'anti- momentum'. 
We have understood the state in terms of two kinds of charges (and their anticharges): 
windings of the elementary string (NSl) and momentum (P). We will call such states 
'2-charge' states, and have found that the entropy of 2-charge states grows as S ^ E. 
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If the string has only left excitations but no right excitations then we get an extremal 
NSl-P state. The entropy is given by setting rip = in (12 ■4p . so we have 5* = 2 y^Tr^ni rip. 
We can use dualities to map this system to other forms. For example we can get D0-D4 
- a bound states of no = ni DO branes and n4 = Up D4 branes. A further T-duality along 
a direction in the D4 gives Dl — D3, where the rii Dl branes are perpendicular to the 
ns D3 branes. We depict this in FigJS] Note that each Dl brane gets 'broken up' into 
ns pieces. Thus there are nin^ fractional Dl branes, and their different positions give 
~ uiUs moduli in a classical description of the branes. Quantizing the wavefunctions on 
this moduli space will again give the 2-charge entropy S = 2\f2Tx ^fn^nl. 
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Figure 3: Different kinds of branes 'fractionate' each other, giving a large entropy. 



2.2 Three charges 

Can we get an entropy that grows with energy faster than S ~ El Let us recall the 
microscopic description of 3-charge black holes. 

Consider type IIB string theory, and let there be 5 compact directions, which we write 

as X S^: 

Mg,! ^ M4,i X X ^1 (2.9) 

We will wrap branes on the compact directions, and obtain an object that is a black 
hole in the in 4+1 nonconmpact directions. The black hole in [7] was made with charges 
D1-D5-P, but since we started with the elementary string above we dualize this to get 
NS1-NS5-P. The NSl branes are wrapped on , the NS5 branes wrap x S*^, and the 
momentum P runs along . The entropy is [7j 

S = 2Ti^nin^np (2-10) 

Let the mass of a brane of type i be m,. Then the energy of the extremal system is just 
given by adding the masses of the branes 

E = niiTLi + nsms + npTUp (2-11) 
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Since the energy of the system is hnear in the numbers of branes n^, we find 

S (2.12) 



Thus the 3-charge entropy grows faster with energy than the 2-charge entropy i\2.7\i . 

Suppose ni,n5 3> Up. Let us add a small amount of extra energy without changing 
the charges, so that the system is no longer extremal. One finds that the Bekenstein 
entropy of the near-extremal hole S{E, Qi) can be reproduced by assuming that we have 
both momentum and anti- momentum excitations, and that these do not interact [8]: 

S = Itx ^ n\n^{^ffip + (2-13) 

E = n^m^ + riimi + (n^ + Upjirip (2-14) 

If only one charge is large, ^ ^i, Up and the system is slightly off extremality, we find 
that we can again reproduce the Bekenstein entropy exactly [15j by writing 

S = 2n^/n^{^yT^ + ^/W^){y/n^ + ^yfl^) (2-15) 

E = n^m^ + (ni + ni)mi + {rip + np)mp (2-16) 

The expression (12.151) needs some explanation. The number is determined by the 
given NS5 charge for the system, but there are four other numbers: ni,ni,np,np. The 
net NSl and P charges [fii^hp) give 

Til — fii = fii, Hp — Up = hp (2-17) 

and a third relation comes from the energy (12.161) . This leaves one free parameter, and 
we should extermize 5* over this parameter to obtain the entropy. The result then tells 
us how the energy wants to partition itself into fractional excitations, and is found to 
exactly reproduce the Bekenstein entropy of the system. 

We can make a natural extension of the above formulae for entropy to the case where 
all charges are comparable, and the system is not close to extremal. This case will 
therefore include the Schwarzschild hole. We write [16] 

E = {n^ + h5)m^ + (ni + ni)mi + {up + hp)mp (2.19) 
Again we have the conditions 

^5 ~ ^5 = 111 — hi = hi, rip — hp = hp (2.20) 

with the energy given by (I2.19P which assumes no interaction energy between the branes 
and antibranes. This time there are 6 parameters and 4 conditions, and we must again 
extremize 5* over the remaining 2-parameter family to get the correct S. The resulting 
S{E, Qi) agrees exactly with the Bekenstein entropy of black holes in 4+1 dimensions, 
for all values of charges Qi and energy E. 

If we scale up all charges and the total energy the same way, we find that for all these 
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entropies arising from using three kinds of charges we get S ^ E^ . 
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2.3 4 charges 

A similar story is found if we compactify an additional direction, so that we have 

Mg,! -> M3,i xT^xS^ X (2.21) 

Now we have 6 directions on which we can wrap objects, and we get black holes in 3+1 
noncompact dimensions. We can take the NS1-NS5-P charges that we had above and 
add a fourth charge: KK monopoles that have 5*^ as their non-trivially fibred circle. We 
can again dualize these four charges to a variety of forms. A form that looks symmetric 
in the four charges consists of four D3 branes that wrap the 6 cycles of the compact 
as follows 

D3u3 D3u, D3246 DS^^e (2.22) 
Thus any pair of D3 branes shares one common direction. The entropy is again given by 

m 

S* = 27r(v/^+ v^)(v^+ v^)(v^ + v^)(v^+ v^) (2-23) 
where we extremize over the nj,ni subject to 

ni-ni = hi (2.24) 

and 

E = J2i^i + ^i)^i (2-25) 

i 

This entropy agrees exactly with the Bekenstein entropy S{E, Qi) of a hole in 3+1 di- 
mensions. 

By changing the orientation of one of the D3 branes we get a nonsupersymmetric 
but still extremal system, and the entropy of this system was matched to the Bekenstein 
entropy recently in [19]. 

Note that the rii grow linearly with E, so the entropy (12.231) grows with energy as 

5-^2 (2.26) 

2.4 Proposal for entropy in the early Universe 

We have seen that if we take a gas of massless particles we get an entropy S ~ E^d^ . 
We can consider excited string states which give S* ~ -E; we have seen that this system 
can be re-interpreted as a 2-charge system where the two charges fractionate each other 
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and produce the entropy. With three charges we get S ~ E2. With four charges we get 
E ~ E^. 

Now consider the Universe, where we have compactified all the spatial directions to 
a torus. The traditional big bang picture envisages a radiation filled Universe at early 
times. In |[2j a string gas was considered and in |3l H] a 'brane gas' was taken. Such gases 
can give entropy S ^ E. But we have seen above that general fractional brane states 
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can give a much higher S at large E. We wish to adopt an equation of state for the early 
Universe that will reflect this high entropy and thus correspond to the most generic state 
for given E. 

We will assume that the entropy has the form 



N 



S = A'l[{^i + y/Wl) (2.27) 



i=l 

Since we will take the net charge of each type to vanish, we have rii = fti and we get 

N N 



S = 2^ A' n = ^ n ^ (2.28) 



i=l 1=1 

The energy is 

E = ^mi{n, + ni) = 2'^mini (2.29) 

i i 

Here mi is the mass of the brane of type i 

m, = Tp\[Lj (2.30) 
j 

where Tp is the tension of a p-brane and the product runs over all the spatial directions 
of the brane. 

We assume that the system is in thermal equilibrium, so we will maximize the entropy 
f l2.27p for given eE To find the state with maximal entropy at a given time t, the Lj 
are held fixed (which fixes the mj), and the total energy is held fixed at E. Taking into 
account this energy constraint we maximize 

N 

S = S- X{Ebranes - ^) = A JJ - A(2 ^ m,n, - E) (2.31) 

i=l i 

Extremizing over rii gives 

= = 2i^ ^'-''^ 

Note that the energy is equipartitioned among all types of branes, each type getting 
energy (there is no sum over k) 

E . , 

Ek = nkiTLk = — (2.33) 



■^To see if the assumption of equilibrium is true, we will have to compute the rate of interactions 
between fractional branes. This interaction depends on the total number of branes in the bound state. 
We do not address these issues here, and hope to return to them elsewhere. 
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2.5 Stress tensor 



We have seen that the entropy of black holes is reproduced by assuming that the energy 
gets partitioned optimally between different kinds branes and antibranes. In this com- 
putation the energy is taken to be just additive; i.e. there was no energy of interaction. 
In [16] it was shown that with this same assumption of noninteraction between branes 
we can reproduce the pressures exerted by the black hole on the various compact cycles. 
Thus on the one hand we can take the black hole geometry and for compact directions 
Ui look at the asymptotic fall-off of gy^yi] this is related to the pressure components of 
the stress tensor in a weak gravity situation. On the other hand we can take the set of 
branes and antibranes that we obtained by extremizing an expression like fl2.18p . fl2.23p . 
compute the pressure each brane exerts by itself on the compact directions, and just add 
these pressures. One again finds exact agreement between the black hole result and the 
microscopic computation]^ We will thus also use a simple sum over the pressures of the 
branes describing our configuration. 

Let us first compute the stress tensor of a single p-brane. The action of the brane is 



S=-Tp J V-^*"'^ rfP+i^ (2.34) 

where g^jf is the metric induced on the worldvolume. The stress tensor is given by 

= (2-35) 

Let the length of the direction be Lj. Let the brane be wrapped on directions . . . x^. 
The volume of the brane is = IliLi The volume of the directions transverse to the 
brane is Vtr = Y[f=p+i ^i- The total volume of the torus is = VpVtr- The stress tensor 
has only diagonal components. We find (there is no sum over k) 

D-l 

T(p)\ = -Tp l[6{x,-x,), k =l,...,p 

i=p+l 

T(^^ = 0, k =p + l,..., (/}-!) (2.36) 

where 6 is the covariant delta function (J S{x)^/—gxx dx = 1), and Xi give the position 
of the p-brane in the transverse coordinates. 

Now suppose there are rip branes of this type, smeared uniformly on the transverse 
directions Xi, i = p + 1 . . . {D — 1) . Then we get 

T(^)\ = -Tp^ = -Tp^ = (2.37) 
Vtr y ^ 



''in [TB] the variables compared between the two computations were certain hnear combinations of 
the pressures; for a direct computation of pressures from wrapped branes see for example (20j . 
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where Ep = UpTpVp is the total energy carried by this type of brane. Using fl2.33p we 
have 

where p = ^ is the energy density. Including the contribution of the corresponding 
antibrane, we get from this type of brane the pressure 

P=-^P (2.39) 

Now suppose there were Ni types of branes wrapping the direction Xj. Then the pressure 
in the direction Xi will be 

P^ = -^P = (2.40) 

where we have defined 

= (2.41) 

Momentum modes P contribute —1 to A^j (they have a positive pressure while branes have 
a negative pressure). Note that the appearing in fl2.4ip counts the types of branes 
wrapping different cycles, not the number of branes along those cycles. An example 
might make this clearer. Take 11-dimensional M theory; then there are 10 compact 
spatial directions. Consider the charges 

M5i2345, M5i2367, M5i4567, A (2.42) 

where the subscripts indicate the directions along which the branes wrap (Pi is momen- 
tum along the direction xi common to all M5 branes). There are 4 kinds of charges; 
thus = 4. Along xi we have the contribution +1 from each of the M5 branes and the 
contribution —1 from P; thus we have A''i = 3 — 1 = 2, and Wi = — ^ = — |. For X2, 
we have a contribution +1 from the first and second types of M5 branes, so we again get 
N2 = 2 and ^2 = — |. A similar result holds for x^, . . . ,xi. Along xg, ■ ■ ■ ,xiq we find no 
charges, so A^j = and the corresponding Wi vanish. Thus we get 

{wi,...,wio} = w = 0,0,0} (2.43) 

2.6 Comparison between our approach and brane gas models 

What new features can strings and branes bring to the early Universe? In [2] the idea 
of 'string gas' was examined; in [3111] the extension to 'brane gases' was considered. We 
will also be using branes, and some of our computations will resemble those in brane gas 
scenarios. But there is a significant difference in the basic idea between our approach 
and brane gases. Here we outline some points of this difference; this discussion should 
help us put forth our conjectured picture more clearly. 
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(a) In a string gas we can wrap strings along arbitrary cycles of the torus. Similarly, 
in the brane gas of [4J M2 branes were wrapped on all cycles {ij) of the torus. But 
in our computation we need to take a set of branes that are mutually BPS; it is only 
such sets that manifest fractionation and the consequent large entropy fl2.27p . M2 branes 
are mutually BPS if they share no common directions, M5 branes need to share three 
common directions, and an M2 - M5 combination must share one common direction. 

(b) In the brane gas model of [1] the entropy comes from vibrations living on the 
surface of the M2 branes. This entropy is proportional to the area of the brane, and at a 
suitably high temperature the energy cost of each unit area of the brane is balanced by 
the entropy carried by that area. This gives a Hagedorn phase with S ^ E. 

In our case the energy E goes to creating certain sets of branes and antibranes that 
fractionate each other, and thereby create the entropy (12.271) . Even if we have just three 

3 

kinds of charges, this entropy grows as S" ~ £"2, much faster than the Hagedorn entropy. 
Thus at high energy densities we expect to get the fractional brane state rather than a 
Hagedorn state. 

(c) In the brane gas model when two branes intersect they tend to annihilate. Thus 
sets of branes that do not generically intersect are expected to last for longer times, and 
govern the long time dynamics of the system. 

The situation is quite the opposite in our case. Consider the three charges NS1-NS5-P 
which we considered for the 3 charge black hole. The NSl is bound to the NS5, and thus 
'lives in the plane of the NS5'. The P charges are carried by excitations of the NS1-NS5 
bound state along the direction common to the NSl and NS5. For fractionation to occur, 
all charges must 'see' each other. 

One may then wonder why the branes and antibranes do not immediately annihilate 
to radiation. But we have already seen that if E is high then there is more entropy in 
the fractional brane state than in radiation, so there should not be an annihilation to 
radiation. Let us discuss the annihilation of brane anti-brane pairs in more detail. 

If we place a brane and an antibrane together, we get a system that can be described 
by a tachyon sitting at the top of its potential hill [21]. Classically the tachyon can sit 
at the top of the hill indefinitely, while quantum mechanically it will fall to its ground 
state. If we take a large number N of branes, and one antibrane, then we can describe 
the branes by their gravity dual. We then find that the antibrane just falls down the 
throat created by the branes, and no radiation emerges for long times [22]; thus there is 
no quick annihilation between the branes and the antibrane. 

In fact the annihilation process for fractional branes has been well studied in the 
black hole context, where one finds that decay of brane-antibrane pairs gives Hawking 
radiation. For the three charge system described by (12.131) we can compute the rate of 
annihilation of PP pairs to radiation, and find that the rate matches Hawking emission 
exactly in spin dependence and grey body factors [P|. Similar agreement is found for the 
4-charge analogue of (12.131) [23] and for the system described by (I2.15P [23]. We can even 
get the exact emission rate for the general hole described by (I2.18l) . (l2.23p . if we boost 
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the neutral hole to add charges; this maps the neutral hole to the near-extremal system 
used in the above mentioned results [25]. While boosting in a compact direction is not 
an exact symmetry of string theory, it may be a good approximation for large charges, 
and is similar to the idea used in Matrix theory. 

All these computations suggest the following picture for black hole microstates. The 
state has a large number of fractional branes and antibranes, and the potential describing 
this system has a large number of saddle points which give metastable states. The system 
slowly drops from one metastable state to a lower energy one, giving Hawking radiation, 
which is a process suppressed by powers of h for classical sized black holes. We thus expect 
that our system on the torus will be composed of fractional branes and antibranes 
with the high entropy fl2.27p . and branes and antibranes will not annihilateEl 

In summary, let us take an analogy from nuclear physics. At low energy we see 
hadrons, but at high density and pressure we get a quark-gluon plasma, where deconfine- 
ment has liberated the elementary degrees of freedom to generate the highest possible 
entropy. At very high energies these elementary constituents are essentially noninteract- 
ing quanta. In our case we have a high energy density in the early Universe, and black 
hole physics suggests that the most entropically favored configuration is one of fractional 
branes. Black hole computations also suggest that these fractional brane quanta are free 
to leading order, and that we should find the total energy and pressure by adding the 
contributions from each brane in the state. 

3 Einstein's equations 

We take the metric to have the form 



The coordinates Xi are compactified with period unity (0 < < 1). The nonvanishing 
components of the connection are 



D-l 




(3.1) 



i=l 




ti ~ 




The relevant components of the Einstein tensor are 




(3.3) 




'''Brane-antibrane models for black holes were also considered in |26| . 



12 



= ^ + ^(y^)_^_y^ + G\ (3.4) 

(There is no sum over k in (13 ■4p .) The Einstein equations are G^i, = SttGT^^. The 
nonvanishing components of the stress tensor are 

T\ = -p, T'^k =Pk = Wkp (3.5) 

so we get the field equations 

i i ' 

ak^ak^ya^^_al_ya^^ g^^^^ _^ ^^^^ ^g^^ 
ak ak ^ ai at ^ ai 

Substituting (13. 6p in (13.71) we get 



+?(E^)-§-e!!^=(i+-')i|(E!^)'-^e|i (3- 



ak _|_ ak 
ak ak '—^ a 



4 A Kasner type power law solution 

For the empty Universe with toroidal compactification we have the Kasner solutions |27j . 
where the radii grow as powers of t. A power law solution was also found for the case of 
isotropically wrapped branes in [3]. We will see that with the equation of state that we 
have chosen we can get a power law solution for any choice of the Wi which characterize 
the brane wrappings. 
Thus write 

ai = tti t^' (4.1) 

Thus 

_ ^ dj _ f3i{(3i - 1) 
Substituting in (13. 8p gives 

l(E.A)-il 



We write 



Then we have 



(4.2) 



^--[ ^ (/-I) ^'-'^ 
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Let us define 

W = Y,w^, U = J2w'^ (4.6) 

i i 

We can get two consistency conditions from (14.51) . First we sum over k in (14.51) . getting 

1R_l1/|2_/1 lR_l42 

J2fl.^A^(D-l)l^^i^^^]-Wl^^^] (4.7) 

Next we square tlie (3^ and tlien add: 

(4.8) 

One solution to tfiese equations is A = 1, B = 1, wliicfi gives tlie well known vacuum 
Kasner solutions [27]. To find other solutions, note that eq. (14.71) is linear in B, and gives 

S . A ^'^-^j,"_f!-/-^' (4.9) 

Substituting in (14.81) we get a quadratic equation for A. Solving this, we get two additional 
solutions, one of which is A = 0. Collecting all these solutions we have the following cases: 



A = 0, B = (4.10) 
This gives /Si = for all i, and thus corresponds to empty Minkowski space. 

(ii) 

A = l, B = l (4.11) 

These are the known vacuum Kasner solutions. Thus there will be no matter, and 
the different expansions of the different directions give a self-consistent solution of the 
Einstein equations. All j3i satisfying (14. 4p with A = B = 1 give allowed solutions. 

(iii) 

^ 2iD-l-W) 



{D~1) + {D-2)U -W^ 
^ (D - 1) + (D - 2YU -2W - {D - 3)W' 

[{D-1) + {D-2)U -W^]^ ^ ■ ^ 

This gives a solution with a nontrivial stress tensor contributed by branes. From (14.51) 
we find 

r 2(iy-l) 2{D-2) 

W - D{U + 1) + 2U + V ^W^ - D{U + l) + 2U + V ^ ' 
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5 The equations in the general case 

Let us write 



Thus 



Equation (13.81) gives 



a,; 



ii + ii 



Let us define 



V 
Q 



i 



li 



li 



(5.1) 
(5.2) 



(5.3) 



(5.4) 



Then (EJ]) is 
Summing (15.51) over k gives 



|7. + 7.P-|P-Q = >^-Q)(1 + -.) 



{D - 2)jV + -{D-1)Q = ^{V^ - Q){D -1 + W) 



Multiplying ( 15. 5p by 7^ and then summing over k gives 

Multiplying (15.51) by Wk and then summing over k gives 

j^S + rS- WjV -WQ = ^{V^ -Q){W + U) 



(5.5) 



(5.6) 



(5.7) 



(5. 



Interestingly, we find that even though there are D — 1 variables 7j, the three moments 
(15. 4p form a closed system of three first order equations. We can write (I5.6p - (I5.8I) in a 
more convenient form by defining 



Q = Q-V^ 



(5.9) 



15 



Then our three equations become 

V + = -KiQ (5.10) 

^ + PQ = -SQ (5.11) 

S + VS = K2Q (5.12) 

where 

- 2{D-2) ^^■^^> 
1 1-W 

K2 = -2^^^W + U] (5.14) 

If V, Q,S are known then we get the 7j from 05.51) 

1 ~ 1-W 

The ttj are then determined by (15. ip . 

5.1 A more convenient form of the equations 

The left hand sides of (I5.10p - (l5.12p have a similar form, which suggests that we define 
an integrating factor. Consider eq. (l5.12p . We can write it as 

4(e4^'^*5)= 7^26^4 (5.16) 
at 

where to is an arbitrary constant that we will take as the initial time where we specify 
data. For any quantity F we write 



Then (15.160 becomes 
Similarly, eq. (l5.10p becomes 
Eq. (l5.1ip becomes 



F = e^4^'^*F (5.17) 
jS = K2h (5.18) 



iv^-K^Q (5.19) 



^fi-sh (5.20) 

Thus in this equation there appears the quantity 5* and not S. Our goal is to get a closed 
system of equations in the hatted variables. To this end we note that for the number 
unity we can write the hatted symbol 



f _ C Vdt H r' Vdt ir 01 ^ 

1 = e-'^o ■ 1 = e-^*o (5.21J 
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Using j we can write fl5.20p as 

|s^-f6 (5.22, 

Note that 

±UV (5.23) 

SO we finally do have a closed system of equations in the hatted variables. We collect 
these equations together for later use 

= -K,Q (5.24) 

at 

d ^ S ^ 

Q = Q (5.25) 

dt I 

= K2Q (5.26) 

at 

It' = * <^-^'' 

Eq. (l5.15p for the 7i can also be written simply in hatted variables 



j^lk = -SkQ (5.28) 

where 

1 1-W 

^'=2^^^ + "^'^ (5.29) 

5.2 Integrals of motion 

The hatted version of the basic equations allow us to note some simple integrals of the 
equations. 

From (15.241) and (I5.26P we find immediately that 



dt 

which gives 



^ {K2V + KiS) = (5.30) 



where the constant is determined by initial conditions. 
From fl^:^ and we find 



K2 - 

S = ——— V + constant (5.31) 



*^' = -^'2 = ^5^P (5.32) 
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which gives 

% = ^V + Fu (5.33) 
where are constants determined by initial conditions. Note that 

E| = 1 (5-34) 

k ^ 

Since '^klk = see that we must have 

^Ffc = (5.35) 

k 

5.3 Physical ranges for parameters 

Note that 

j i 

where V is the volume of the spatial torus. We have to choose a direction of time to call 
positive, and we can use this freedom to require that at the time to where we give initial 
conditions 

V{t = to) > (5.37) 
The integrating factors that converts un-hatted quantities to hatted ones is 

/ = e4^'^* = ^ (5.38) 

Thus throughout the physical range of evolution we will have 

/ > (5.39) 

From fl5.17p it follows that hatted and un-hatted variables have the same sign. 
From (13. 6p we find that 

^(E^^)' -\Y.^" = - 2) = ^^GP (5-40) 

i i 

SO the energy density is 

Our matter is made up of the quanta in string theory, and we have seen that the energy 
of different kinds of objects will be simply added to obtain the total energy. Each of 
these quanta have a positive energy, so we will have 

p>0, Q<0, Q<0 (5.42) 
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The total energy is E = pV, so Q is just the total energy upto a (negative) constant 

S = ^S = -^J5 (6.43) 

We have found that pk = WkP- We will assume the dominant energy condition, which 
states that for each direction the pressure satisfies pk < p. (As in the case of energy 
density, the pressure is obtained by adding the pressure contributed by the different 
string theory objects in our state, and these satisfy the dominant energy condition.) 
Thus we have 

Wi < 1 (5.44) 
for all i. This gives W < (-D — 1), or equivalently, D — 1 — W > 0. This implies that 

From fl5.42p and (15.451) we find that the RHS of (15.241) is non-negative, so P is a 
non-decreasing function of time 

> (5.46) 

dt - ^ ^ 

From (I5.37P and (I5.46P we see that for all t > to we will have 

V>0 (5.47) 

From (I5.27P we see that 

so that / = ^ is a nondecreasing function of time. Thus the Universe will not 'recoUapse' 
to 1/ ^ 0. 

The variables S, S can have either sign, and this sign can change during the evolution. 



6 Solving the equations 

We observe that in the system (I5.24p - (l5.27p . three of equations have Q on the right hand 
side. We can divide by Q and absorb it in the definition of time, by writing 

-^^^4- (6-1) 
dt dr 

We have put in the negative sign because Q is negative; with this sign, the variable r 
increases when t increases. Thus the t, t variables are related by 

T=fdt'{-h), {t-to)= r ^ (6.2) 
J to Jo (_g) 
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where we have chosen the lower hmit of t to be the time to where we specify initial 
conditions. The system fl5.24l) - fl5.27l) then gives 

IV ^ (6.3) 

= i (6.4) 

= -K, (6.5) 

dr 

d ^ V , , 

— / = -— (6.6) 



dr 

We can immediately solve (16.31) and (16. 5p : 



Q 



V = KiT + Ai (6.7) 
S = -K2T + A2 (6.8) 



where Ai,A2 are constants. Now (16.41) and (16. 6p become 



d X S {-K2T + A2) , , 

— Q = ^ = ^ 6.9 

dr I I 

A| . (6.10) 
dr Q Q 



From these equations we deduce that 



:-^Q)i + Qi^i) = ^{Qi) = -{K, + K2)t + {A2 - A,) (6.11) 
dr dr dr 



which gives 



QI = -{K, + K2)- + {A2 - A,)t + A, (6.12) 



where ^3 is another constant. Taking / from this equation in substituting it in (16.91) 
gives 

= t^3l±A^ t (6.13) 

dr -(Ki + K2)^ + {A2 ~ Ai)t + As 

or 

:^[log(-&] = t^^l+A^ (6.14) 

dr^ ~{K^ + K2y-^ + {A2-A^)T + As ^ ' 

where we have written (— Q) in the argument of the log since Q is negative. The quadratic 
in the denominator on the right hand side can be written as 

- {K, + K2)y + {A2 - A,)t + A, = JJ^1±M (^ _ r^){T - r^) (6.15) 
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where ri, r2 are the two roots of the quadratic. We now need to know if these roots are 
real or complex, and if they are real, then where r lies on the real axis with respect to 
these roots. We will study these roots below, but for now we write the formal solution 

to (Km 

(— Q) = y44(r — ri) {^fi+-'f2){'-i-'-2) (r — (■^i+-^2)(n-'-2) (6.16) 

where is a constant. 

Thus all variables V, Q, S, I have been expressed algebraically in terms of the time 
parameter r. From fl6.7p we see that upto a suitable choice of origin and a constant 
scaling, r is just the variable V. Thus if we use V to measure time, then all other variables 
are given by rational functions of this time. To get back to the physical problem however, 
we need to relate r to t. This is done through (16. 2p 

1 r , 2(~riK2+A2) 2(-r2K2+A2) 

(t-to) = — / (r' - ri) (^fi+^2)('-i-'-2) (r' - r2) WT+K^WT^dr' (6.17) 
^4 Jo 

The integral on the RHS is given by an incomplete Beta function. This function is defined 
by |28] 

B,{p,q)= I s^'\l-sy-Hs (6.18) 
Jo 

The precise expression for (16.171) in terms of the incomplete Beta function will depend 
on the location of r with respect to the roots ri, r2. 

Since the relation between t and r is transcendental, we will analyze the solutions 
qualitatively to see the dynamical behavior that results for different choices of parameters. 

6.1 Different dynamical behaviors 

Consider the integral in (I6.17p . For what follows we recall eq. (15.460 which says that V 
cannot decrease with time. There are three possible cases: 

(a) The integral (I6.17P diverges at a finite value of r. Then we reach t = 00 with 
finite r. Then from (16. 7p we see that V asymptotes to a finite constant. 

(b) The integral (16.170 diverges as r ^ 00. In this case "P ^ 00 as t ^ 00. 

(c) The integral (16.170 converges. In this case we have a divergence P ^ 00 at a 
finite time t. 



6.2 Dependence on parameters 

We now wish to see which of the above behaviors results for which choices of parameters 
and initial conditions. Recall that / is positive (eq. (l5.39p ) and Q is negative (eq. (l5.42p ). 
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Thus the left hand side of (16.121) is negative. Consider the function 

fir) ^ i-IQ) = i^^l±Mr^ _ - A)r - A, (6.19) 

The physical values of parameters then requires 

/ > (6.20) 




(a) (b) 

Figure 4: (a) Downward facing parabola for Ki + K2 < (b) Upward facing parabola 
for Ki + K2 > 0. In each case a physical choice of parameters leads to motion along the 
bold line segment. 

6.2.1 Ki + K2<0 

In this case the parabola is concave downwards. From (16.201) we see that only the part 
of / above the r-axis describes a physical evolution. If the roots ri,r2 in (16.151) are real, 
then the parabola will intersect the r-axis, and a part of the parabola will lie above this 
axis (Fig|l](a)). If the roots are complex, the parabola will be entirely below the r-axis. 

In this case it is easy to show that ri,r2 will be real. The discriminant of the poly- 
nomial in (I6.15P is 

A = {A2 - A,y + 2A3{K, + K2) (6.21) 

^Here and in other computations below we consider only generic values of the parameters for sim- 
plicity. For example we do not explicitly look at the border Ki + K2 — between the two cases below; 
such special cases can be easily worked out explicitly. 
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Note that ^3 is the value of the negative quantity QI at the initial time r = 0, so 

A3 < (6.22) 

Since Ki + K2 < in the present case, we find that all terms in A are positive and thus 
A > 0. Thus the roots ri,r2 are real, and the parabola will look as in FiglD^a). The 
evolution will take place on the solid part of the parabola in FigHl^a). Since the evolution 
ends at a finite value of r, we find from (16. 7p that V asymptotes to a constant, and thus 
we will be in case (a). 

The initial data is given at r = 0, which must be a point between the two roots of 
the parabola. Thus ri < and r2 > 0. Since ri < t < r2 during the physical evolution 
the relation (I6.17P should be written as 



it-to) = r^ I (r' - nr (^2 - rTdr' (6.23) 







where we have defined 



2i-r,K2 + A2) 2{-r2K2 + A2) 

(K^ + K2){r^ - r^y ' {K, + K^^n - r^) ^ 



We find 



(t-to) = ^(^2 -ri)^i+^2 (Br-r, [ai +1,^2 + 1] - + 1,02 + 1]) 

(6.25) 

In Fig|5] we plot graphs for an example that illustrates case (a). We let Wi = 
{.9, — .9, — .9, — .9, — .9, — .1, — .1, — .1, — .1, — .1}. We have set to = 2 and have taken 

liit = to) = 1, ai{t = to) = 1 for S'll ^- We plot V,—Q, and one from each set 

having the same Wi. Note that — Q is proportional to the total energy in the Universe. 



6.2.2 Ki + K2>0 

In this case the parabola is concave upwards. With a little more effort we can again show 
that for physically allowed initial conditions, the discriminant A is positive, and thus 
ri,r2 are real; this computation is done in Appendix flA.ip . Thus the parabola intersects 
the r-axis, as shown in Fig|l](b). From (16.201) we see that physical motion can take place 
only on the two segments of the parabola above the r-axis. Thus there appear to be two 
possible branches, a 'left' branch and a 'right' branch. With some effort we can show 
that for a physical choice of parameters, we cannot be on the left branch; this is done in 
Appendix (lA.2p . Thus motion takes place only on the right branch, as indicated by the 
solid part of the parabola. 

We see that the motion will extend to r — > 00; thus from (16.70 we see that we will 
be in case (b) or in case (c). To distinguish between these cases consider the integral 
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fl6.17p . Since we are on the right branch above the r-axis we see that (r — ri) > and 
(r — > 0. Thus from (16.161) we find A4 to be a real positive constant. In (I6.17P we 
find that the large r behavior gives 



2Ko 



t-to~ j dT'{T') (r) K,+K2+' (6.26) 

Thus the convergence of this integral is determined by the sign of /i = —j^^^ij^ + 1- 
Note that we have Ki + K2 > in the present part of the analysis. So we can equally 
well ask for the sign of {Ki + K2)fi = Ki — K2. We thus have the two cases 

(i) Ki — K2 > 0: In this case the integral (16.171) diverges at r = 00, and we have 
case (b). 

(ii) K1 — K2 < 0: In this case the integral (16.171) converges and we have case (c). We 
note however in Appendix (]A.3l) that we can have this case only if at least one of the Wi 
is less than —1. It is conventionally assumed that Wi lie in the range —l<Wi<l. The 
upper limit comes from the dominant energy condition, but there is no strong reason to 
require the lower limit. For the quanta that we get from string theory though we do have 
— 1 < Wi < 1, as can be seen from the definition (12.411) . 

In either of these cases (b),(c) we are on the 'right branch' of the parabola in FigJD^b). 
Thus the point r = where the initial data is specified lies to right of the two roots of 
the parabola. So ri < r2 < 0, and we find from (I6.17p 



- rij°i+"2+i / s 
(t-to) = — — fir-ra [0:2 + 1, -tti - «2 - 1)] - -Bi:! [^2 + 1, -«! - a2- 1)] 

I/I4I V T-ri ri J 

(6.27) 

In Figl6] we plot graphs for an example that illustrates case (b). We have taken 
Wi = —.2 for all i. We have set to = 2 and have taken 7j(t = to) = 1, ai{t = to) = 1 for 

all i. We plot V, — Q, and ai. 



6.3 Solving for 7^, 

From (I5.32P we find 

which gives 

where fk are constants. Since 



at dr 



Q 



% = SkT + fk (6.29) 



5^7fc = /$^7fc = /^ = ^ (6.30) 
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Figure 5: Plots of V,—Q and a selection of for Wi — 

{.9, -.9, -.9, -.9,-. 9,-. 1,-.!,-. -.1}, a set that gives + K2 < and 
illustrates case (a) behavior. We have taken 7^(2) — 0^(2) — 1 for all i. We see that V 
asymptotes to a constant. 
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Figure 6: Plots of V, Q,ai for the choice Wi = —.2 for all i. This gives Ki + K2 > 0, 
Ki > K2, and thus case (b) behavior. V grows without bound. 
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we have one relation between the fk 

5^/fc = ^(r = 0)=Ai (6.31) 

k 

Now note that 

— (logafc) = — = = ^ = (6.32) 

at a-k I I 

Thus 

— (logafc) = (--) -(logafc) = — = -[ -] (6.33) 

where we have used fl6.12l) . fl6.15p . This gives 

Qk =Ck {r - ri) iK:,+K2Kn-r2) (^T- _ ^2) WP^imn^ (6.34) 
where Ck are constants. 



6.4 Evolution as a function of t 

We have seen above that all variables can be expressed as algebraic functions in terms of 
the time parameter r. But in the metric (13. II) the natural parameter is t, and the relation 
between r and t was given through the incomplete Beta function, so this relation is not 
easy to picture. We can however find that evolution as a function of t at late times, 
where the relation between r and t simplifies. 



6.4.1 Case (a) at large t 

From FigJU^a) we see that the evolution takes us towards the root r = r2. In Appendix 
(IA.4p we show that 

aa + 1 < (6.35) 
Then we see from (16.231) that t diverges as r — r2 

t ^ (r2 - r)°^+^ (6.36) 

and from (I6.34p we get 

2(gfcr2+/fc) 

Ofc ~ i2(-r2A-2+A2)-(ifi+K2){ri-r2) (6.37) 

Thus the late time behavior of the depends on the initial conditions that specify 
parameters like fk and Ai,A2,A3 which give the roots ri,r2. In Appendix (1A.4I) we 
observe that the power of t in (I6.37P can be written in a simpler form 

ttk ~ t T(-2) (6.38) 
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6.4.2 Case (b) at large t 

In this case we have t — > oo as r — oo. Let us examine this large t behavior. From 
(ESSD we get 

r~t^i-^2 (6.39) 

From fl6.34p we find 

ttk ~ ~ t^*^ (6.40) 

where we have used the expression f l5.29p for 6k and the expressions (I5.13p . (l5.14l) for 
Ki, K2 to recognize that the power of t is the same as the power appearing in (14.131) 
for the Kasner type solution. Thus we see that evolution for case (b) asymptotes at late 
times to the Kasner type solution for the given Wfclll 

6.4.3 Case (c) for t tf 

In this case the r integral converges 

2(-r-^K2+A2) 2(-r2K2+A2) 

(r - ri) (^fi+^f2)('-i-'-2) (r - r2) (^i+^2){'-i-'-2) = - (6.41) 



\M Jo 

To find the behavior as t ^ we note from (16.1 7p that in this limit r is large and we 
have 

T {tf -t)^i-''2 (6.42) 

From fl6.34p we get 

2(5^ 

ak (tf -t)T^T^ {tf -tf'' (6.43) 

where we have noted that the powers of {tf — t) appearing here are the same as those 
that appear in the Kasner type power law solution ~ t^. 



7 Discussion 

We have assumed an equation of state suggested by black holes; the energy goes to 
creating 'fractional branes' which have a high entropy and can thus dominate over other 
kinds of matter. If we assume that the state of the Universe is always the maximal 
entropy state for the given energy E, then we get an equation of state Pi = Wip, for 
which the dynamic can be solved analytically. 

We discussed some of the ideas behind the fractional brane state in section (12. 6p . If 
we have a few branes in a given volume of space, then these branes will tend to annihilate 
to massless quanta. But if we increase the energy E in the given volume to very large 
values, then it becomes entropically favorable to produce a large number of suitable sets 

^In [4j it was also observed (from numerical computations) that there was a kind of 'attractor mech- 
anism', so that solutions with generic initial data became similar to each other at late times. 
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of branes and anti-branes. Branes in such a set are mutually BPS, and 'fractionate' 
each other, producing an entropy that grows more rapidly with E than the energy of 
radiation or a Hagedorn type string or brane gas. As we discussed in section ( \2.6\} . black 
hole results suggest that at these densities the fractional brane quanta seem to behave 
as essentially free objects. Thus the energy E and pressures pi for the state are given 
by just adding the contributions from the branes and antibranes. The fractional branes 
do not find it easy to find each other and annihilate; the rate of this annihilation can be 
computed for black holes where it reproduces exactly the rate of Hawking radiation. 

If we follow the state of the present Universe backwards in time, the energy density 
keeps increasing. It has been postulated that for sufficiently early times we reach a 
Hagedorn phase of strings. This phase has no pressure, but does have energy; thus as 
we go further back in time E does not change {dE = —PdV) but a 7^ 0, so the density 
grows still further. This suggests that at sufficiently early times the fractional brane state 
should be present. 

We have not analyzed the astrophysical implications of the evolution we find; there 
are many issues to be addressed and we hope to carry out a detailed study of these 
elsewhere. Here we outline some of these issues and raise some relevant questions. 

We have not taken any specific choice of the Wi] rather we solved the problem for an 
arbitrary set of Wi. Which choice of Wi that gives the largest possible entropy for M theory 
on T^°? We can find large sets of mutually BPS branes with 10 compact directions, but 
we need to prove that choosing branes and antibranes of these varieties will indeed give 
the entropy fl2.27p . We thus need to generalize the brane constructions that have worked 
so well for black holes and understand the entropy of high density states in string theory. 

For generic choices of Wi we get a power law expansion, and for many choices the 
power is too low to give us some kind of 'infiation'. Note than in infiation we find the 
infiaton in a low entropy state, while in string/brane gas approach we seek the maximal 
entropy state for the given energy. In this sense we are closer to the latter approach; we 
look for a maximal entropy state but observe that the Hagedorn gas is not the highest 
entropy state in string theory. 

Infiation gives a reason for the sky to look homogenous on large length scales. But 
there is a different source of long distance correlations that can arise with the fractional 
brane gas. In black holes semiclassical analysis suggests that quantum gravity effects 
are confined to planck or string length, but because of fractionation we find that the 
brane bound state has a size that grows with the number of quanta in the state, and 
thus we get nonlocal quantum effects all across the interior of the horizon [22]. So in 
our Cosmological fractional brane state we can also expect quantum correlations across 
macroscopic distances. 

What is the fate of the fractional branes as the Universe expands? Consider the 3- 
charge entropy fl2.18p and the 4-charge entropy fl2.23p . For a given energy E, is it more 
advantageous to create three kinds of charges or four? We see that if all the exceed 
unity, then the 4-charge entropy is higher, but if one of them drops below unity, then 
the 3-charge entropy will be higher. Such a transition was studied in [TTl I2U] where it 
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gave a microscopic description of the black hole - black string transition. In our present 
problem the mass of each type of brane changes as the evolve, and it is possible that 
some Hi drops below unity. In that case we would have to continue the evolution with a 
different set of branes and thus a different set of Wj. 

Clearly it is crucial to know how many branes form the fractional brane state; this 
determines the rij. Thus if the Universe was infinite and all branes were linked up to 
form the fractional brane state, then each would be infinite and the entropy per unit 
volume would diverge. (Note that doubling the volume more than doubles the entropy if 
three or more types of charges are involved.) It is clearly entropically advantageous for 
more and more of the matter to be linked up in the fractional brane state, but when the 
Universe starts the islands of matter so linked are presumably small; they might then 
grow rapidly as different islands come into causal contact. 

Note that the density of matter in a fractional brane state need not be high as long as 
enough total matter is linked up in the state. For example we can make a black hole with 
arbitrarily low density matter as long as there is enough total energy E in the ball of 
matter. Any fractional brane matter left over today could show up as a dark component 
with its own dynamics. 

All these are interesting questions, and we regard the present work as just a first pass 
on the problem with these ideas; we hope to return to the above issues elsewhere. 
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A Conditions arising from physical ranges of param- 
eters 

First we collect together some definitions and relations that will be used to establish the 
inequalities in later sections of the Appendix. 

Our spacetime dimension is D] thus the number of space directions is 



d = D-l 



(A.l) 



We can regard the Wi as a vector in d dimensional space 



(A.2) 



Similarly 



7 = {7i,72,---,7d} 



(A.3) 
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It will be convenient to define another c?- dimensional vector which has all components 
unity 

r={i,i,...,i} (A.4) 

Then we find 

W = l-w, U = w-w (A. 5) 

and 

P = 7-l, 2 = 7-7, S = 'J ■ w (A. 6) 

Note that 

\w\ = Vu, \l\ = Vd (A.7) 

Let the angle between 1 and w be 6. 

Let the initial value of 7 be 70. The vectors 1 and w define a plane in d dimensional 
space. We can decompose 70 into a part along this plane and a part perpendicular to 
this plane 

7o = 7|| + 7± (A.8) 

Let the angle between 7|| and 1 be 0. 

Note that the equations fl5.24l) - fl5.27p are invariant under the rescaling 

V^fJ^V, Q^fi^Q, S^fiS, t^^-H (A.9) 

where /i is a constant. This rescaling just corresponds to scaling time, which scales the 
dj while keeping fixed the a^. Thus the initial 7^ can be scaled without changing the 
essential behavior of the evolution. We will use this freedom to set 

|7||| = 1 (A.IO) 

This will simplify our expressions without loss of generality. We have 

(l-7o) = (1-711) = v^cos0, (w.70) = (w.711) = v^cos(6'-0), (7o-7o) = 1 + |7± 



With a little algebra we find from the definitions fl5.13l) . fl5.14p 

Ki + K, = --^^-—^{[{D-1)U-W']-[U + {D-1)-2W]) 

^ — {[Ud -W^]~[U + d- 2W]) 



2 

(A.ll) 



2{d 

^^^^ ([(1 ■ r)(^ • w) - il-w)'] - [(1 - 1) + {w-w)- 2{l.w)]) 

(A.12) 

^— (^[t/c/sin^^] - [t/ + d-2v^cos^]) (A. 13) 



2{d 

Using fl6.4p . fl6.6p we get for the function / giving the parabola fl6.19p 

4zm = -(i4-Q+Q4-l]='P-S (A.14) 



dr \ dr dr 
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A.l Proof that A > for the case Ki + K2 > 

Recall that 70 = 7|| + 7±. We can write 

%=al + Pw (A.15) 

At the initial time t = towe see from (15.171) that hatted variables equal the corresponding 
un-hatted variables, and from fl5.38p we note that /(to) = 1- Thus using (lA.llll 

(70-1) = P(r = 0) = P(r = 0) = (A.16) 
i^-w) = 5(r = 0) =5(r = 0) = ^2 (A.17) 

Taking the inner product of (1A.15I) with 1, w we get two equations, solving which we get 

A1U-A2W „ -A,W + A2d 

Note further that 

(70 ■ 7o) = Q(r = 0) = (Q + V^)iT = 0) = (S + V^)ir = 0) = A, + Al (A.19) 
We now compute 

If-i-T = \io-(yl-Pw\^ 

Al(U(d-l) -W^) - Ald + 2AiA2W , 

Since |7±P > 0, we get 

^ Al{U{d-l)-W')-Ald + 2A,A2W 
^ - ^ Ud-W^ ^ ^ ■ ^ 

Note that in the present case we have Ki + K2 > 0, so we can multiply the above 
inequality by Ki + K2 without reversing the sign of the inequality. We thus have 

A = {A,-A2y + 2A,{K, + K2) 

> M o..^ , ,,, Al{U{d-l)-W')-Ald + 2A,A2W 

> iA,-A2) -2{K, + K2) jjj-^, 

[A, {U[d-l) + W{l-W)) + A2{W-d)f 

{d-l){Ud-W^) ^ ■ ' 

Note that 

{Ud-W^) = {w ■w)il-l) - il-wf > (A.23) 
by the Schwartz inequality. Thus we have 

A > (A.24) 
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A. 2 Absence of left branch of parabola for Ki + K2 > 

In Fig|4](b) positivity of energy density requires that we be above the r-axis. Above this 
axis we find a left branch and a right branch of the parabola, and we wish to show that 
the left branch is not allowed for a physical choice of parameters. 

We will use four constraints that together will rule out the left branch: 

(i) In FigJ7|(a) we have drawn the vector 1 as OA, and the vector w as OB. The 
angle between these vectors is 6. Let L be the length of BA. Let L_i_ be the length of the 
perpendicular from O to the line BA. 

Note that in the present case we have Ki + K2 > 0. Using ( 1A.12I) this gives 

[{l-l){w-w) - {1-wY] < [{1 ■ 1) + {w ■ w) - 2(1 ■ w)] (A.25) 
Let the area of the triangle BOA be Aa- The left side of is ( ]A.25[) is 

(1- l){w-w) - {1-wy = (VuVdsine)^ = {2A^y = {LL^y (A.26) 
The right hand side of flA.25|) is |(1 — w)P = L^- Thus flA.25|) gives us the constraint 



< 1 (A.27) 

(ii) From the dominant energy condition we have Wi < 1 for all i. This gives 
W = Y.i'Wi < d, so 

d-W = {l-w)-l>0 (A.28) 

(iii) From fl5.4ip and the positivity of the energy density p we have at r = 

Q = Q-V^ <0 (A.29) 
From flA.6|) and using flA.lOl) we get 

l + l7±r-^^(r = 0) <0 (A.30) 

which gives 

V{t = 0) > 1 (A.31) 

(iv) Since we are looking at the left branch of the parabola given by the function 
/(r), we have ^/ < 0. Using (1A.14P we find 

S>P (A.32) 
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Figure 7: Graphical representation of physical constraints for Ki + K2 > 0, left branch 
of parabola. 



center O; the tip of 7|| must lie on this circle. We draw the tangent AB' to this circle, 
touching the circle at C. 

(i') The lines AO and AB' form a wedge. There are two parts of this wedge: one 
to the left of A and one to the right; the latter is shaded. From simple geometry using 
the condition (i) (eq. (lA.27|) ) we find that the point B must lie inside this wedge. 

(ii') The condition (ii) tells us that B cannot lie in the shaded part of the wedge. 
This follows because the line BA is the vector 1 — w and OA is the vector 1. For B in 
the shaded part of the wedge the angle between BA and OA is obtuse, and thus (1A.28P 
is violated. 

(iii') Recall that the angle between 7|| and 1 is (p. The condition flA.3ip gives 



Now we will use these four constraints. Since 




V = (1-70) = (1-711) = Vdcoscf) > 1 



(A.33) 



This tells us that 



\(t>\ < (t>m = COS 



-1 



1 



(A.34) 



Thus the tip of 7|| must lie in the arc CC. 



(iv') The condition (EH is 



(1 — tu) ■ 7 = (1 — w) ■ 7|| < 



(A.35) 
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This tells us that the tip of 7|| must lie in the arc DD'D" . 

Finally, simple geometry tells us that if B is indeed in the wedge indicated then the 
arcs CC and DD'D" do not overlap. This proves that the four constraints (i)-(iv) are 
incompatible, and the left branch of the parabola is not obtained for physical values of 
parameters. 



A. 3 Impossibility of case (c) for —l<Wi<l 
From the Schwartz inequality we have 

{l-l){w -w) - {1-wf >0 (A.36) 

which gives 

d 



U - — >0 (A.37) 



Using ( 1A.37I) we find 



A little algebra gives 

U-^ = 2{-K^ -K, + 2^ K^) (A.38) 

K2 < -Ki + 2^-^ Kf (A.39) 
a 

Recall from fl5.45p that Ki > 0; this followed from the dominant energy condition Pi/p = 
Wi < 1 which we have assumed. Thus we can divide by Ki in the above inequality 
without reversing the inequality signs 

^ < -1 + 2—— K, (A.40) 

Now suppose that we have 

Wi>-1, i,l,...d (A.41) 

Then 

W = J2^i>-d (A.42) 

i 

Then from the expression (15.13^ for Ki we find 

D-l-W d-W d 

^^ = ^(d3^ = 2(Z3I)<Z^ (^-43) 

Using this in flA.40p we have 

< 1 (A.44) 
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Noting again from (15.451) that Ki > 0, we can multiply through in this inequality without 
reversing the inequality sign 

Ki-K2>0 (A.45) 

Since case (c) arises for Ki — K2 < 0, we see that to get this case we will need to violate 
(1A.41I) for some Wi. 

Note that if we have Wi = —1 for all i then we get a Cosmological constant term in 
the Einstein equations, and thus we get exponential evolution. Though this is a faster 
expansion than the power laws we found for case (b), we still do not get a singularity at 
finite t. 




A. 4 Absence of singularity at finite t in case (a) 

The physical part of the graph in FigJD^a) ends at a finite value of r, but we would like 
to know if this corresponds to a finite or infinite value of the physical time parameter t. 
We will show that for physical values of parameters the point t = T2 cannot be reached 
at finite t. 

First we show this by a simple argument using the equation (15.251) . which gives 

t c 

-^dt (A.46) 

to I 

Recall that I{t = to) > 0, and from (15.271) and (15.471) we know that it cannot decrease. 
Thus the denominator of j is bounded below. On the other hand from (16. Sp we see that 
\S\ is bounded for our interval ri < r < r2. Thus j is bounded for our evolution from 

r = to r = r2. But from (16.121) we see that at the root r2 we must have Q = (since / 
is bounded below and cannot vanish). So we see that the left hand side of (1A.46P must 
diverge for r = r2, but from the right hand side of (IA.46|) we find that this divergence 
cannot occur at any finite t. 

Now let us establish this behavior directly by proving (I6.35p . which is 

2(-r2K2 + A2) 



[K^ + K2){n-r2] 



+ 1 < (A.47) 



In the present case we have Ki + K2 < 0, and from our ordering of roots we have 
Ti — r2 < 0, so we can multiply the inequality (IA.47P by {Ki + K2){ri — r2) without 
reversing the inequality sign 

- 2{-r2K2 + A2) + {Ki + 7^2) (ri - r2) < (A.48) 

We write 

V^[r = ri)=Vu Vo{r = r2)=V2, So{r = n) = Si, So{r = = S2 (A.49) 
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From (I621),(I6SD we get 

- r2K2 + A2 = S2, K,{n - r2) = Vi - A, K2{r, - r2) = S2 - Si (A.50) 

and the inequality to be established becomes 

(A - A) - {Si + 4) < (A.51) 

Consider the parabola /(r) drawn in FigJD^a). At the two roots ri,r2 the slope of 
/(r) must be equal and opposite. Thus using flA.14p we have 

iVi-Si) = -{V2-S2) (A.52) 

Using this in flA.Sip we find that the inequality to be established is 

- 2P2 < (A.53) 

We have already seen in fl5.47p that V is nonnegative throughout the evolution. We will 

get V{t2) = only if Q ~ p = everywhere, which corresponds to empty Minkowski 
space and thus gives no singularity anywhere. For 7^2 > we get (1A.53I) . which establishes 
(1X471) . 

Using the notation above we can obtain fl6.38p . From fl6.29p . fjA.50|) . flA.52l) we see that 



2(4'"2 + fk) 27fc 7fc 7fc 



2{-r2K2 + A2)~iKi + K2)iri-r2) 2S2 - {Vi - V2 + S2 - Si) A ^2 

(A.54) 

which gives f l6.38p . 
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